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Abstract We review X-ray constraints on dark matter in giant elliptical galaxies 
(10 12 M@ < M vir < 1O 13 M ) obtained using the current generation of X-ray satel- 
lites, beginning with an overview of the physics of the hot interstellar medium 
and mass modeling methodology. Dark matter is now firmly established in many 
galaxies, with inferred NFW concentration parameters somewhat larger than the 
mean theoretical relation. X-ray observations confirm that the total mass profile 
(baryons+DM) is close to isothermal (M ~ r), and new evidence suggests a more 
general power-law relation for the slope of the total mass profile that varies with the 
stellar half-light radius. We also discuss constraints on the baryon fraction, super- 
massive black holes, and axial ratio of the dark matter halo. Finally, we review 
constraints on non-thermal gas motions and discuss the accuracy of the hydrostatic 
equilibrium approximation in elliptical galaxies. 



1 Introduction 

Cold, non-baryonic dark matter (CDM) is a critical ingredient of the widely ac- 
cepted Big Bang cosmological paradigm (QHHIil). Although the dark matter (DM) 
particle has yet to be directly detected in a terrestrial laboratory El 0), the total 
amount of DM is known precisely from observations of the large-scale geometry 
of the Universe, in particular the cosmic microwave background radiation, distant 
Type la supernovas, baryon acoustic oscillations, galaxy cluster mass functions and 
baryon fractions. These observations indicate that DM comprises about 23% of the 
energy density budget of the Universe today compared to about 4% for baryonic 
matter QUI)- 

While the standard cosmological model describes the large-scale Universe very 
well, it has encountered interesting tension with observations of DM on galaxy 
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scales. Perhaps most notably, beginning with the seminal studies of rotation curves 
of dwarf spiral galaxies dHO, it is now generally recognized that both small and 
large disk galaxies have central radial DM profiles that are flatter than initially pre- 
dicted by the standard model when considering only the gravitational influence of 
the DM {Ilklll El. 

Similar cored DM profiles are obtaine d by op tical studies of 
elliptical galaxies at the center of galaxy clusters ( 14 ; 15 ; 3 17 ; 18 ), though X-ray 
observations of some clusters find "cuspy" total mass profiles consistent with the- 
ory This "core-cusp" problem has generated great interest as it may yield 
information on the nature of the DM that is apparently inaccessible to observations 
on larger scales. An influential early suggestion to reconcile observation and theory 
is that the DM particles are self-interacting (12 It |22t 23t 24h. th ough observations 



20 



apparently do not currently favor that possibility (Il9t 

Alternatively, the core-cusp problem may reveal clues to the influence of baryons 
on galaxy formation which should be most important near the centers of galaxies. 
In the standard model slightly over-dense regions in a nearly uniform early DM 
distribution acted as seeds for future galaxies and clusters. Baryonic material sub- 
sequently fell into the potential wells established by the DM. The interplay of the 
baryons and DM through, e.g., adiabatic contraction and dynamical friction, has 
likely altered the spatial distributions of both quantities during galaxy formation 
and evolution (H3 M, EH E3 E ; I32II33I). Hence, the DM profiles on galaxy scales 
can serve as important laboratories for studying galaxy formation and the nature of 
the DM particle. 

Like all the different types of galaxies, giant elliptical galaxies have their own 
special properties that make them deserving of dedicated study. Unlike disk galaxies 
or the smallest dwarf elliptical galaxies, giant elliptical galaxies are sufficiently mas- 
sive that multiple powerful techniques to measure their mass profiles are available 
for many systems. In particular, only for giant ellipticals can gravitational lensing 
and X-rays be used to complement the information available from stellar dynam- 
ics, the latter of which is also available for smaller galaxies. Also, in contrast to the 
smallest galaxies, every giant elliptical is thought to host a super-massive black hole 
(SMBH) at its center, the formation of which is apparently intertwined with the host 
galaxy d34t 1351) . It is also becoming increasingly recognized that galaxies acquire 
their gas in two modes that are distinguished by whether or not the in-falling gas 
is heated (3^; 37)- Both modes can operate in the transition mass regime of giant 
elliptical galaxies (10 12 - 1O 13 M ) $3$). 

One of the advantages of using X-ray observations to probe DM and SMBHs in 
elliptical galaxies is that the hot X-ray-emitting interstellar medium (ISM) fills the 
three-dimensional galactic potential well, thus providing a continuous tracer from 
the galactic nucleus out to well past the optical half-light radius where stellar dy- 
namical studies become very challenging. In order to translate an X-ray observation 
of the hot ISM into a gravitating mass profile and DM measurement, it is required 
that hydrostatic equilibrium holds to a good approximation throughout the region 
of interest; i.e., that the thermal gas pressure of the ISM balances the weight of the 
gas. With X-ray observations providing increasing evidence of morphologically dis- 
turbed elliptical galaxies, it is timely to review the accuracy of the approximation of 
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hydrostatic equilibrium in X-ray determinations of their DM profiles. Understand- 
ing and quantifying deviations from hydrostatic equilibrium not only will lead to 
more accurate X-ray mass determinations, but any measured non-hydrostatic gas 
motions will provide clues to elliptical galaxy formation d39l l36l [371: [38 : 4(J 41). 

Since giant elliptical galaxies tend to be located in dense environments, it can be 
difficult to disentangle any DM associated with the elliptical galaxy from its parent 
halo. Naturally, most of the attention of X-ray studies has been devoted to those 
systems with the largest X-ray fluxes which has lead to an emphasis on elliptical 
galaxies located at the centers of the most massive galaxy groups and clusters. While 
all elliptical galaxies are of interest for DM studies, in this review we address galaxy- 
sized halos by focusing on the lower mass regime that is still accessible to detailed 



X-ray study (10 M < M v ; r < 10 M ) We restrict consideration to objects where 
a single giant elliptical galaxy dominates the stellar light, which at the upper end of 
the mass range under consideration are often classified as fos sil grou ps d42t l43t k4) 
which are thought to be highly evolved, relaxed systems (145tl46l : |47|). 

In the present review we focus on the constraints on DM in elliptical galaxies 
obtained using the current generation of X-ray observatories - Chandra, XMM, and 
Suzaku; constraints from the previous generation of X-ray satellites are reviewed by, 
e.g., (I48tl49l : l39l) . While we take care to mention relevant results obtained from other 
techniques, our review is not intended to be complete from the multi-wavelength 
perspective. The interested reader is urged to consult the reviews by ( 5(J 51 ) and ( 52 ; 



53l) for discussions of the constraints on DM in elliptical galaxies obtained by stellar 



dynamics and gravitational lensing. 



2 Hydrostatic Models of the Hot ISM 

2.1 Preliminaries 

The hot ISM of an elliptical galaxy (T ~ 10 7 K - 0.9 keN/k B ) has properties that 
are extremely similar to the (hotter) intracluster medium (ICM) of a massive galaxy 



cluster (154t l55t 1561) . While the detailed radiation spectrum is described well by a 



collisionally ionized plasma of electrons and ions obeying the coronal approxima- 



tion (1571) . the macro-structure of the hot ISM in an elliptical galaxy is for many 
purposes also described well by a simple monatomic ideal gas. We highlight four 
key features of the "hot gas" relevant for constructing hydrostatic models to measure 
the DM content in elliptical galaxies. 



Collisional Fluid 



In a fully ionized plasma, the equipartition timescale for electrons considering only 
their mutual Coulomb interactions is (1581) . 
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where m e is the electron mass, e is the electric charge, kg is Boltzmann's constant, n e 
is the electron number density, and T e is the electron temperature. For T e > 4 x 10 5 K 



the Coulomb logarithm of an electron-proton plasma is (15 81) . 
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The mean free path of the electrons for Coulomb collisions is v e x e - e Sis El, 
where v e = \fik^f e ~Jm e is the rms thermal velocity of the electrons, so that, 
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Because m e does not enter the above equation this expression also represents the 
mean free path for the protons (for n e = n p and T e = T p ). The electron and pro- 
ton mean free paths are generally much smaller than the length scales of interest 
for DM studies indicating that the plasma may be treated as a collisional fluid. 
Still, in the outer regions near the virial radius, where n e is small, and therefore 
the Coulomb mean free path is large, it is possible the fluid approximation con- 
sidering only Coulomb collisions may not always be justified. But the hot ISM of 
elliptical galaxies may typically possess tangled ~ 1/iG magnetic fields J62ll63l). If 
so, the electron gyroradius, r g = v e /o, for angular gyrofrequency (0 = eB/m e c, 

— ,z-.(^)" 2 ( S )" 2 (^)- pc , „ 

is vastly smaller than the relevant length scales for DM studies. (The same is true 
for the proton gyroradius obtained by setting m = m p .) Consequently, if the hot 
ISM contains (plausible) weak magnetic fields, then it is expected to behave as a 
collisional fluid on all length scales of interest for DM studies, allowing an isotropic 
pressure and temperature to be defined for any fluid element. 



Local Maxwellian Velocity Distribution 

In any fluid element in the hot ISM Coulomb collisions between electrons establish a 
Maxwellian velocity distribution with kinetic temperature T e on the time scale given 
by equation (fl}. This is much shorter than other time scales associated with the hot 
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ISM in a relaxed elliptical galaxy. In general the ions can reach equilibrium with a 
kinetic temperature different from that of the electrons. But electron-ion collisions 
will establish equipartition on a time scale ss (m p /m e )T e - e « 1870T e - e « 10 7 yr for 
the conditions typical of the hot ISM d58l) . Only for galaxies currently undergoing 
strong dynamical disturbances, particularly in the central regions of some galax- 
ies with AGN, might be evolving dynamically on a comparable time scale so that 
the electron and ion temperatures could be different. In relaxed elliptical galaxies, 
however, the hot ISM is expected to behave locally as a Maxwellian velocity distri- 
bution with a single kinetic temperature for the electrons and the ions. Therefore, the 
thermal pressure of the hot ISM is characterized by that of a monatomic ideal gas, 
PgasksT I jJ,m a , where is the mean atomic weight of the gas and m a is the atomic 
mass unit. For a fully ionized gas of solar abundances, n = 0.62. 



Thermal Pressure Dominates Magnetic Pressure 

Observations of radio halos and relics in several massive galaxy clusters indicate 
that the ICM is magnetized with a field strength of 1 - lOjliG (1641) . Unfortunately, 
elliptical galaxies are not observed to possess the large-scale diffuse radio emission 
observed in some massive clusters, and thus knowledge of their magnetic fields is 



derived from systems with embedded radio jets where the inferred ~ fiG fields d62l) 
may not be representative of the typical hot ISM. For fields of this magnitude, with 
^mag = B 1 j%% and P theim = n gas k B T « l.9n e k B T, 

-0.015 fJLWiV'f - ,5) 
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which implies a minor, though very uncertain, contribution of magnetic fields to 
the pressure support of the hot ISM. Many theoretical models of inflowing gas in 
elliptical galaxies predict that small seed magnetic fields are amplified possibly to 
this level or larger d39l and references therein), although the small gas fractions 



observed at the centers of groups (I65t 1661) likely indicates that any hypothetical con- 
centration magnetic fields, e.g., arising from the amplification owing to inflowing 
gas just described, has been dispersed by feedback. X-ray studies of DM in ellip- 
tical galaxies generally neglect the magnetic pressure, an assumption that can be 
tested indirectly by (1) comparing observed stellar mass-to-light ratios to theoreti- 
cal stellar population synthesis models and (2) comparing mass profiles obtained by 
X-ray methods with those obtained from stellar dynamics and gravitational lensing 



Thermal Motions Dominate the Hot ISM in Relaxed Galaxies 



When the gravitational potential evolves on a time scale longer than the dynamical 
time, approximately given by the crossing time of the stars, 
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where d is the length scale of the region under consideration and a is the stellar 
velocity dispersion within that region, a time-independent description of the grav- 
itational potential is appropriate. This quasi-static description of the gravitational 
potential is valid regardless of the dynamical state of the hot ISM, because the 
gas contributes negligibly to the total mass (gas+stars+DM) over most of the ra- 
dial range within the virial radius of an elliptical galaxy. Radiative cooling of the 
hot gas acting alone cannot drive evolution on time scales smaller than the cool- 
ing time («5x 10 6 /« e yr for T = 10 7 K and solar abundances, or 5 x 10 9 yr for 
n e = 10~ 3 cm~ 3 ) which is generally much longer than f cr0 ss- But rapid evolution on 
a time scale < f cr0 ss is expected to prevail in the central regions of galaxies that have 
experienced recent, strong AGN feedback presumably generated in response to gas 
cooling. For our purposes, a relaxed elliptical galaxy is one in which the structure 
of the hot ISM is also evolving quasi-statically on a time scale longer than f cross 
described by a state of hydrostatic equilibrium, which applies provided any non- 
thermal gas motions (i.e., motions that are not associated with the random velocities 
that define the thermal gas pressure) are negligible compared to the sound speed, 



1/2 f T xl/2 

where 7=5/3 is the adiabatic index. This implies a sound crossing time similar 
to the stellar crossing time in equation (0. Turbulence is expected to be introduced 
into the hot ISM from merging and intermittent AGN feedback. Interesting limits on 
turbulent motions have been obtained via constraints on resonance scattering in the 



center (< 1 kpc) of the elliptical galaxy NGC 4636 (1671) and from Doppler broaden- 
ing of the X-ray emission lines in a small number of elliptical galaxies (NGC 533, 
NGC 1399, NGC 4261, NGC 5044) from the sample of Q suggesting that any 
turbulent motions present in the central regions of those systems are highly sub- 
sonic. While not expected to dominate, turbulent velocities may not be negligible 
since cosmological hydrodynamical simulations typically find that turbulence con- 
tributes up to ~ 20% of the pressure support of the ICM in relaxed galaxy clus- 
ters (69; 7(| 71; 72; 73). The irregular, subsonic Ha velocities observed at the cen- 



ters of some elliptical galaxies also suggest turbulent velocities in the hot gas (74 
It is unlikely that rotation dominates the ISM dynamics on large scales in elliptical 
galaxies given the lack of flattened disks observed (T75h . though modest rotational 
spin-up might occur in the innermost (< 1 kpc) regions of some galaxies q4 II). The 
expected inflow velocities arising from radiative cooling are negligible (1391) . and a 
recent study of galaxy-scale halos using a cosmological hydrod yna mical simulation 



concludes that their hot ISM is generally quasi-hydrostatic (I76tl77l 
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2.2 The Equation of Hydrostatic Equilibrium: Usage Guidelines 

For the relaxed elliptical galaxies under consideration, hydrostatic equilibrium of a 
fluid element of hot ISM is expressed as the balance between forces per unit volume 
of gravity and gas pressure, 

VPgas - -pgasVO, (8) 

where <P is the gravitational potential, P gas is the pressure and p gas the density of 
the ISM. If non-thermal effects are significant that can be represented by a pressure 
(P nt ), such as for isotropic random turbulence or a magnetic field, then P gas = P t +Pnu 
where P t is the thermal pressure. For the case of solid-body rotation, <P can be re- 
placed by an effective potential. For general rotation (and other ordered gas motions) 
the convective derivative term (v ■ V)v in the Euler equation must be added to the 
left-hand side of equation ||8}. As remarked previously, the hydrostatic approxima- 
tion is appropriate provided the non-thermal gas motions are subsonic. 

In order to measure directly the magnitude of any non-thermal motions present in 
the hot ISM, it is necessary to measure precisely the Doppler shifts and broadening 
of emission lines which is beyond the capability of the current X-ray satellites - ex- 
cept possibly for the innermost regions of some galaxies noted above. Consequently, 
indirect methods are required to assess whether a galaxy is suitably relaxed for the 
hydrostatic approximation. Assessment of the morphology of the X-ray image is a 
convenient first step to locate candidate relaxed galaxies for the following reason. 
Because the total underlying stellar and DM distribution is expected to have a shape 
close to that of an ellipsoid, the gravitational potential is also expected to be ap- 
proximately ellipsoidal, though rounder than the generating mass distribution. Since 
in hydrostatic equilibrium the X-ray emissivity traces the same three-dimensional 
shape as the potential for any temperature profile also see q6b, galaxies 

that exhibit irregular, non-ellipsoidal features in their X-ray images might not be 
suitably relaxed for hydrostatic study. The existence of such asymmetrical features 
(e.g., substructure, cold fronts) does not guarantee large departures from hydrostatic 
equilibrium, as demonstrated by cosmological hydrodynamical simulations of clus- 
( 69t 70t 71 : 72 : 73). especially when the features are sufficiently localized so 



ters 



that they can be easily excluded from study 



Another important consideration in hydrostatic analysis is to insure that the X-ray 
emission that is attributed to hot ISM is not, in fact, heavily contaminated by unre- 
solved discrete sources. A large number of the brightest discrete sources (mostly 
low-mass X-ray binaries - LMXBs) can be detected and excluded in nearby galax- 



ies using the high spatial resolution of Chandra (see|81j, and references therein); the 
lower resolution of XMM leads to far fewer detections, and only a few sources (if 
any) can be detected with Suzaku. Because the combined spectrum of unresolved 
LMXBs is spectrally harder than the hot ISM, the spectral signature of the hot ISM 
can be reliably extracted from the data provided it is not overwhelmed by the unre- 
solved source component. The situation is more problematic in X-ray faint galaxies 
(having low ratios of X-ray-to-optical luminosity) like NGC 3379, where hydro- 
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static mass estimates have been attempted (182c 1831) but the soft, diffuse X-ray emis- 
sion attributed to hot ISM most likely originates almost entirely from cataclysmic 
variable stars and stellar coronae (84; 85]). 

Once a suitably relaxed galaxy with sufficient emission from hot ISM has been 
identified, the accuracy of the hydrostatic equilibrium approximation can be as- 
sessed indirectly via consistency checks. First, one may examine if any well- 
motivated hydrostatic model is able to provide an acceptable fit to the observed 
density and temperature profiles of the hot ISM. If an acceptable fit is obtained, the 
data are consistent with, though do not necessarily require, hydrostatic equilibrium. 
A poor fit - particularly one with sharp discontinuities in the density and tempera- 
ture profiles - provides strong evidence for significant violations of the hydrostatic 
equilibrium approximation^. Second, the gravitating mass profile inferred from the 
X-ray data can be compared to independent measurements using other techniques 
having different assumptions (stellar dynamics, gravitational lensing). When results 
obtained from different methods agree, it provides support for the underlying as- 
sumptions of each method and implies that systematic errors are well understood 
(see $£}. 

Hence, given the limitations of current X-ray data, a sensible approach to apply 
and assess hydrostatic equilibrium is the following: 



Guidelines for Applying Hydrostatic Equilibrium in Elliptical Galaxies: 

• Select systems which tend to have a regular, approximately cicular or elliptical 
X-ray image morphology. If the image has asymmetric features, they should be 
preferably confined to spatial regions that are small compared to the region of 
interest; e.g., the inner region showing evidence of possible AGN feedback. 

• Unresolved discrete sources should not dominate the X-ray emission. As a rough 
guide, at least half the X-ray luminosity within the optical half-light radius should 
originate from hot gas. 

• Determine whether a hydrostatic model is able to provide an acceptable fit to the 
density and temperature profiles of the hot ISM over the region of interest. 

• When possible, compare the mass profile obtained via the hydrostatic approx- 
imation with that obtained by an independent method; e.g., stellar dynamics, 
gravitational lensing. The hydrostatic equilibrium approximation is judged to be 
useful if it provides DM measurements of comparable (or better) quality than 
other techniques - or provides any information when no other technique is avail- 
able. 



The disturbed X-ray image of M 84 is a spectacular example ( 89) of the profound 



effect that AGN feedback can have on the hot ISM and is not a galaxy where hydro- 



1 Still, even in such cases the hydrostatic approximation might be useful. The core of M 87 - a 
system more massive than we consider in this review - possesses various spatial and spectral (I86t 
[871) irregularities that are very likely associated with sizable non-hydrostatic gas motions. Yet 0881) 
found that the gas is close to hydrostatic away from the regions of most significant disturbance. 



Dark Matter in Elliptical Galaxies 



9 




static equilibrium is expected to apply very accurately. But not all elliptical galax- 
ies display such strong, large-scale irregularities in their X-ray images. In Figure Q] 
we display several elliptical galaxies with mostly regular images suitable for hydro- 
static analysis. Two of the them, NGC 720 and NGC 4649, have very regular images 
deep into their cores (Figure |2). Another two, NGC 4472 and NGC 4261, possess 
low-level image fluctuations consistent with cavities associated with radio jets that 
do not greatly interfere with hydrostatic mass analysis. For several of these galax- 
ies, radial entropy profiles have been computed (19 0t l91t 1921) and found to increase 
monotonically with radius as also found in groups and clusters (93), indicating that 
the hot ISM is stable to convection and consistent with approximate hydrostatic 
equilibrium. 





Fig. 2 Residual significance images of the central regions of NGC 720 (top, adapted from (92b) 
and NGC 4649 (bottom), indicating deviations from a smooth model fit to the X-ray isophotes i90c 
l92l) . Both galaxies appear to be relaxed and symmetrical at all scales accessible to the Chandra 
images. 



2.3 On the Incidence of Relaxed Elliptical Galaxies 



In a controversial study, Diehl & Statler (95h investigated the X-ray and opti- 
cal isophotal ellipticities in a heterogeneous sample of early-type galaxies to test 
whether they are consistent with hot ISM obeying hydrostatic equilibrium. The ex- 
pected correlation between the ellipticities assumes that both hydrostatic equilib- 
rium holds exactly and mass follows (optical) light, so that the stars define the shape 
of the gravitational potential and thus the shape of the X-ray emission. In an attempt 
to mitigate the impact of any DM halo, Diehl & Statler restricted their analysis of 
each galaxy to within a small radius - 60%-90% of the optical half-light radius. 
When the expected correlation between the X-ray and optical isophotal ellipticities 
was not found, they concluded that hydrostatic equilibrium is not ubiquitous and, 
therefore, should never be assumed (even approximately) in mass analysis. They 
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speculate that assuming hydrostatic equilibrium will always result in an error typ- 
ically as large as the value of the mass measured; i.e., relaxed elliptical galaxies 
do not exist. While we agree (as detailed in the previous section) that care should 
be exercised in the application of hydrostatic equilibrium to the study of elliptical 
galaxies, we strongly disagree with their general conclusions about the accuracy of 
X-ray mass measurements. 

The observed lack of a correlation does not allow conclusions to be drawn with 
certainty for an individual object; it only indicates that hydrostatic equilibrium is un- 
likely to hold perfectly in a galaxy randomly drawn from their sample, provided it 
is chosen without any consideration of its morphology. As noted above, hydrostatic 
equilibrium as expressed by equation (|8) is an approximation and is not expected to 
hold exactly given especially the impact of turbulence on the hot ISM induced by 
past mergers and AGN feedback. In other words, no elliptical galaxy is expected to 
be perfectly relaxed. The objects displayed in Figures Q] and [2] are among the most 
relaxed elliptical galaxies possessing high-quality X-ray data. But Diehl & Statler's 
sample also contained galaxies with pronounced large-scale asymmetries that we 
would not recommend the routine application of hydrostatic mass methods follow- 
ing the guidelines discussed previously (e.g., NGC 4636: (19 61) and M 84: (89j)), and 
since they focused on the smallest scales (where X-ray point-source removal is most 
challenging and where AGN-driven disturbances are most serious), in contrast to the 
larger scales that are most important for DM analysis, the implications of their study 
for DM measurements of morphologically relaxed systems are dubious and provide 
no basis for claiming typical mass errors > 100% on all scales in such systems. 

An example of how conclusions regarding the center of a galaxy need not apply 
on larger scales is provided by NGC 5044. In this system Diehl & Statler measured 
an X-ray ellipticity (0.41 ±0.08) greatly exceeding the optical value (0.07 ±0.01), 
which cannot be reconciled easily with hydrostatic equilibrium (unless the gas is 
also rapidly rotating). But the X-ray ellipticity measurement is confined to the re- 
gion within the optical half-light radius (< 4 kpc) containing obvious disturbances 



presumed to be associated with AGN feedback (1971 ; |98j l99l) . (We note, however. 



that there is recent evidence that turbulent motions of the hot ISM are negligible 
compared to the sound speed in the center of NGC 5044 (1681) .) The X-ray emission 
can be traced o ut at least to several hundred kpc, about one-third of the virial ra- 
dius ( lOOl : 101 ). Outside of the very central region discussed by Diehl & Statler, the 
X-ray isophotes are largely symmetric and nearly circular, and acceptable fits with 
hydrostatic models are obtained (1651) . 

Now turning our attention to the central regions within the optical half-light 
radius relevant to the study of SMBHs (fQ, we remark that additional factors 
can introduce scatter into the correlation between the isophotal ellipticities gen- 
erated by the stars and hot ISM. As noted above, Diehl & Statler assume that 
the mass traces the same shape as the optical light, whereas both X-ray and 
optical observations suggest DM distributed differently from the stars typically 
contributes 20%-30% within the half-light radius in giant elliptical galaxies (see 
SjUi. Furthermore, their novel technique to create X-ray images of diffuse hot 
ISM strictly requires there are no spatial variations in the spectral properties; 
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i.e., every galaxy is assumed to be isothermal with constant metallicity through- 
out, i n conflict with obse rva tions of many galaxies possessing high-quality X-ray 
data ([lol[lM[IS[IS[il[IS[Io3[S- Finally, ISM rotation might be impor- 
tant in some galaxies and yet have little relation to stellar rotation because slowly 
rotating inflowing gas at large radius might be spun up to dynamically relevant ve- 
locities if angular momentum is conserved. But rotational spin-up is not guaranteed 
for every galaxy, and instead the gas can be outflowing gently in a subsonic wind 
greatly reducing the impact of any rotation initially provided by the stars - see dis- 
cussion in (1391) . 

Before leaving this topic, we caution that a statistical appraisal of the proper- 
ties of the hot ISM in any sample of elliptical galaxies should consider the sample 
selection criteria. Our impression of the X-ray properties of elliptical galaxies is 
obtained within the context of what are essentially optically selected samples. But 
for the study of ISM properties, it is important to select the sample according to the 
ISM luminosity and not the total X-ray luminosity which also contains the emission 
from discrete sources associated with the stellar population. Diehl & Statler follow 
the common approach of studying a "heterogeneous" sample of elliptical galaxies, 
which essentially amounts to collecting all objects in the archives with data of suf- 
ficient quality. This procedure resulted in a sample containing some well known 
objects with extended, X-ray luminous ISM (e.g., NGC 4472) and some where vir- 
tually all the X-ray emission can be attributed to discrete sources (e.g., NGC 3379 
-HQ). 

It is instructive, therefore, to consider a sub-sample extracted from their sam- 
ple of 54 galaxies that is selected based on the X-ray luminosity of the diffuse gas. 
Using the tabulated diffuse gas luminosities published by Diehl & Statler, we con- 
struct such a sub-sample by excluding massive group/cluster systems (gas L x > 3 x 
10 42 erg s _1 ) and gas-poor galaxies (2(7 lower limit on gas L x < 1 x 10 40 erg s _1 ) 
resulting in a sample of 26 galaxies, give or take a few depending on how we inter- 
pret the quoted lower or upper limits for some systems. This sub-sample includes 
the seven galaxies shown in Figure Q] as well as several other systems with gener- 
ally regular X-ray images and where hydrostatic equilibrium previously has been 
profitably employed - NGC 533, NGC 1399, NGC 1404, NGC 3923, NGC 4555, 
NGC 5044. The sample also contains well-known disturbed systems such as M 84, 
NGC 4636, and NGC 5846, which are expected to deviate from hydrostatic equi- 
librium. In sum, an ISM-selected sub-sample culled from Diehl & Statler's het- 
erogeneous sample indicates that approximately 50% of the elliptical galaxies are 
sufficiently relaxed for hydrostatic analysis. 



2.4 Methods 



The variety of techniques that have been developed to apply hydrostatic equilibrium 
to the study of the mass distributions in elliptical galaxies (as well as galaxy clus- 
ters) can be divided into two broad classes that will be denoted here as smoothed 
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inversion and forward-fitting. Whereas the forward-fitting approach assumes param- 
eterized models for the mass profile and one thermodynamical variable of the gas 
(e.g., temperature profile), smoothed inversion instead introduces one or more ar- 
bitrary parameters to guarantee a smooth (physical) mass profile. The decision of 
whether or not to adopt a particular approach depends on several factors which we 
discuss below. Generally, the smoothed inversion approach is preferred if one is 
primarily interested in measuring the value of the enclosed mass within some ra- 
dius without making any reference to an assumed (parameterized) mass model. If 
instead one desires to test the viability of an assumed mass model (e.g., power-law, 
NFW) and to measure its parameters for comparison to other galaxies and theory, the 
forward-fitting approach is preferred. Whenever possible both approaches should be 
employed to assess the magnitude of any biases. 

We note that in the following it is generally assumed that the gas density and tem- 
perature profiles (projected or deprojected) of the hot gas have been already obtained 
via an analysis of the spatially resolved X-ray spectra. The hydrostatic models then 
use these measured profiles to constrain the mass distribution. This decoupling of 
the spectral fitting and mass modeling, though convenient and generally less com- 
putationally expensive, is not necessary. The formalism of both of these approaches 
is reviewed in Appendix B of (1651) . 

In the following subsections, we outline in more detail how the two broad model 
classes are implemented. 

2.4.1 Smoothed Inversion 
Pros: 

• Mass measured without reference to an input parameterized model. 

• Does not require an input parameterized model for the gas temperature, density, 
pressure, or entropy. 

Cons: 

• Deprojection requires a high degree of symmetry (usually spherical). 

• Cannot self-consistently account for the projection of galaxy emission from any 
radii outside of the observed data range or, generally, extrapolate the hydrostatic 
model outside the data range. 

• Arbitrary smoothing required to obtain a physical (monotonically increasing) 
mass profile. The amount of smoothing increases for lower data quality. 

The aim of the smoothed inversion approach is to measure the mass profile with- 
out reference to an assumed (parameterized) model for the mass distribution or any 
quantity associated with the hot ISM (gas temperature, density, pressure, or en- 
tropy). For this to be achieved, the three-dimensional distribution of the hot ISM 
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must be obtained by non-parametrically deprqjecting the data on the sky. This is 
only possible if the ISM has a high degree of symmetry, and therefore spherically 
symmetric deprojection using the well-known "onion pe el" method originally de- 
veloped for clusters is almost always adopted (1109t ll 101) . Since the onion peeling 



starts from the outer edge of the galaxy and works its way inward, such deprojection 
formally requires that the data span the entire size of the galaxy. As this is not gen- 
erally the case, ad hoc prescriptions are usually employed to estimate t he projection 
of gas from radii larger than the extent of the available data 4ml El. 

Because no global mass model is assumed, the value of the mass determined 
within a given shell can be negative due to noise in the data (amplified by depro- 
jection). To circumvent this problem so as to achieve a mass profile that increases 
monotonically with radius, the inversion procedure must smooth the gas quantities 
(e.g., temperature and density) and the mass profile itself. The amount and type of 
smoothing is arbitrary, and typically as the data quality decreases the amount of 
smoothing required to achieve interesting constraints generally increases. As dis- 
cussed below, types of smoothing include ( 1 ) the use of ad hoc smooth parameter- 
ized models for the gas density and temperature, (2) assuming the ISM is constant in 
spherical shells so that the widths of the shells set the smoothing scales, and (3) one 
or more parameters in a Bayesian prior. While the smoothing can be well-motivated, 
it necessarily introduces additional assumptions, usually expressed by additional 
parameters that must be determined empirically by the data. Consequently, we be- 
lieve that the terms "non-parametric" or "model-independent" sometimes used to 
describe smoothed inversion methods are misleading and should be avoided. 

In the context of DM studies, smoothed inversion is best suited to providing a 
measurement of the enclosed mass within a particular radius without assuming a 
specific parameterized mass model; e.g., measuring M500 for the purpose of estab- 
lishing the mass function of elliptical galaxies. But if one also desires to fit and 
measure the parameters of a mass model (e.g., global power-law slope, NFW con- 
centration parameter), the ad hoc smoothing adds a barrier between the data and 
model which is not present in the forward-fitting approach that is more naturally 
suited for this problem. 



Traditional Approach: Parameterized Density and Temperature Profiles 



For a spherically symmetric galaxy where the hot ISM is supported by ideal gas 
pressure, equation dS) can be solved for the total gravitating mass (stars, gas, and 
DM) as a function of the gas density (p gas ) and temperature (T), 



M(< r) 



rk B T 



jlm a G 



(ilnpga. 

d\ar 



d\nT 
dlnr 



(9) 



where Icq is Boltzmann's constant, jj. is the mean atomic weight of the gas, m a is the 
atomic mass unit , and G is N ewton's gravitation constant. Beginning with M 87 over 
thirty years ago (111 3c ll 141) . this equation has been widely applied to galaxies and 
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clusters traditionally by inserting smooth parameterized functions for the density 
and temperature into equation (0. The ad hoc parameterized functions do not guar- 
antee a physical mass profile that increases monotonically with increasing radius, 
and their use offsets any gain in rigor from not assuming a specific mass model. 
Thus, while this traditional approach has the advantage of being conceptually sim- 
ple and computationally inexpensive, it is best employed as a rough estimate and a 
check on other methods. 



Smoothing by Binning 

The derivatives in equation (0 may instead be computed directly from the density 
and temperature values measured in adjacent radial bins rather than by assuming 
any parameterized models. Because of bin-to-bin statistical fluctuations in the den- 
sity and temperature (or pressure) profiles, increased binning (i.e., smoothing) of 
the radial profile is generally required to produce an approximate monotonically in- 
creasing mass prof ile. This method has only been attempted for the brightest galaxy 



clusters (11 15b ll 161) . since only for data of the highest quality can interesting results 



be obtained without resortin g to heavy binning of the data 



Recently, Nulsen et al. Gil 71) have proposed a scheme that obviates the need 
for explicit computation of the derivatives by assuming that within any spherical 
shell the gravitating matter density and gas temperature are constant. By requiring 
that the gas pressure be continuous across shell boundaries, the gas density is then 
computed within any shell using an isothermal solution of equation ([8]). But the 
pressure continuity condition involves the difference in the gravitational potential 
between adjacent shells, which is related to a derivative (and is the argument of an 
exponential). Consequently, the mass profile derived from this procedure remains 
subject to similar bin-to-bin statistical fluctuations that occur when the derivatives 
in equation (O are computed explicitly. To improve the smoothness of the derived 
mass profile, Nulsen et al. advocate requiring that the gravitating matter density 
decrease monotonically with increasing radius. 



Bayesian Inversion with Smoothing Prior 

A ve ry general approach to smoothed inversion has been propo sed b y Das et al. 
( 118 ), wh ich is based on the seminal works of Merritt & Tremblay dl 19t) and Magor- 



nan Das et al. enforce smoothness via a Bayesian prior which penalizes large 

second derivatives between radial bins in the temperature, circular velocity (hence 
the mass), and pressure gradient. The magnitude of the smoothing is controlled by 
a parameter X, the value of which is arbitrary and depends on both the number of 
radial bins and the model grid spacing. For convenience Das et al. use simulations 
to select a fixed X for all galaxies. Since the sensitivity of the derived mass profile to 
X should be investigated for any galaxy, a reasonable extension would be to define 
a prior for X and then marginalize over X for each system. 
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Pros: 

• No formal restrictions on galaxy geometry (though spherical symmetry usually 
assumed). 

• Allows for self-consistent extrapolation of hydrostatic model and projection of 
galaxy emission to and from any radii outside the data range. 



Cons: 

• Requires an input parameterized model for the mass. 

• Requires an input parameterized model for one of the following thermodynamic 
quantities: gas temperature, density, pressure, or entropy. 

The forward-fitting approach uses parameterized models for the mass and one 
thermodynamic variable of the hot ISM (density, temperature, pressure, or entropy). 
Because the galaxy is completely specified by the models, in principle any geo- 
metrical configuration can be accommodated, though spherical symmetry is usually 
assumed. The model can be extrapolated to radii outside the available data range al- 
lowing for self-consistent treatment of the projection of the ISM of the entire galaxy. 
Forward-fitting can be applied to data that have been deprojected (as with smoothed 
inversion) or directly to the data on the sky. The ability to fit projected models to 
data on the sky is especially important for lower quality, noisy data whose integrity 
would be unacceptably degraded by deprojection noise. 

This approach is preferred when it is desired to test the viability of a particu- 
lar mass model and to measure its parameters. If the assumed model is incorrect, 
then the measured mass at a given radius may differ systematically from the true 
value. It is therefore necessary to explore multiple mass models to gauge the mag- 
nitude of any bias. Since the measured mass may be sensitive to the assumed mass 
model, when possible results using forward-fitting should be compared to those 
from smoothed inversion. 

Prior to Chandra and XMM, when accurate spatially resolved temperature pro- 
files were only available for a few of the most massive galaxies, it was customary 
to make strong assumptions about the temperature profile in order to obtain an in- 
teresting constraint on the DM. Typically the gas was assumed to be isothermal or 
that it obeyed a polytropic equation of state (P gas °< pjas)- Th e most widely used of 



these approaches is the well-known isothermal "beta model" dl21tll22Hl23l) . While 
such procedures still can be useful for obtaining DM estimates from observations of 
low data quality, we focus our discussion on methods intended to exploit accurate 
spatially resolved temperature profiles, as such data are now more widely available. 
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Density-Based 

If an assumption is made about the form of the gas density profile, e.g., that it is 
described by a parameterized function, then taken together with an assumed param- 
eterized mass model the equation of hydrostatic equilibrium (equation [8]l may be 
solved to give the temperature, 



where 7b = T(ro) and p gas ,o = Pgas(/o) evaluated at some reference radius ro. Here 
and below we assume spherical symmetry unless noted otherwise (i.e., in fQ. The 
parameters of the mass and density models and the normalization 7b are then con- 
strained by the observation. As with all the types of forward -fitting approaches we 
discuss, it is common practice first to measure the temperature and density profiles 
from fitting the spectra and then constrain the free parameters of the hydrostatic 
model via a simultaneous fit to the measured profiles. Although the gas contribution 
to the total mass M{< r) in elliptical galaxies can generally be neglected until one 
approaches close to the virial radius, a self-consistent hydrostatic model should in- 
clude it. For these density-based models this is straightforward since the gas density 
profile is assumed, and therefore the gas mass profile is defined before the tempera- 
ture profile is computed. 

An importan t variation of the den sity-based approach was proposed by Fabian 
and colleagues (I109t ll lOt Il24t 1 1251) which exploits the relatively high statistical 
quality of the X-ray surface brightness profile. Sometimes referred to as the "Cam- 
bridge Method", this approach assumes spherical symmetry and analytically de- 
projects the surface brightness using the onion peel procedure to give the three- 
dimensional emissivity profile, e °c Pg as A (T,Z), where A (T,Z) is the plasma emis- 
sivity which depends on temperature T and metal abundances Z and is integrated 
over the relevant (broad) wave band. If A(T,Z) is assumed to be constant over the 
galaxy, then pgasC**) x \/ £ ( r )- Because of the relatively high statistical quality of 
the surface brightness (and emissivity), the density profile derived in this way usu- 
ally can be evaluated in much finer spatial bins than is possible for the temperature 
profile determined from detailed spectral fitting. This allows the finely binned den- 
sity profile to be treated as a continuous function by interpolating values between 
the bins, which can then be inserted (along with a parameterized mass model) into 
equation (ITOb to predict the temperature profile for comparison to the values mea- 
sured from the spectra. 

It is a significant advantage of the Cambridge Method that it does not require for 
input an assumed, possibly oversimplified, parameterized model for the density pro- 
file. Still, the density profile generated assuming a constant A leads to hydrostatic 
solutions that are not self-consistent in the presence of temperature and metallic - 
ity gradients. This effect is not expected to be large because variations in A within 
a galaxy or cluster are probably no more than ~ 20% which translate to inferred 
density differences of only ~ 10%. In principle, the best-fitting temperature pro- 




(10) 
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file could be used to define a new spatially varying A from which a new density 
profile could be computed. The process could be iterated to achieve improved self- 
consistency. 



Temperature-Based 

If an assumption is made about the form of the temperature profile, then in con- 
junction with a parameterized mass profile the equation of hydrostatic equilibrium 
(equation[8]l may be solved for the gas density, 



P gas (r)=pga S ,o^Texp 



jUm a G [>' 1 M{< r) ' 
T(r) r 2 



(11) 



where spherical symmetry is again assumed. This approach can be especially use- 
ful when the temperature profile is poorly constrained by the data. Simple profiles 
(constant, power-law) can be adopted in order to explore their impact on the inferred 
mass profile. Unlike the density-based case, the gas mass contribution to M{< r) 
cannot be included fully self-consistently in the fit. However, because the gas mass 
can be treated as a small perturbation, an iterative procedure can be employed; i.e., 
after first taking the gas mass contribution to M(< r) to be zero, use the resultant 
best-fitting p gas (r) to compute a fixed gas mass profile that is added to M{< r) and 
then fit again. 



Entropy-Based 

Recent studies have suggested that the entropy is the lo gical ther modynamic vari 



able for which to assume a parameterized model (9(1 91 



126U127 ). Consistent with 



the hydrostatic equilibrium approximation is the requirement that the hot gas be sta- 
ble against convection at every radius. Convective stability provides an important 
additional constraint on the mass profile that has be en us ed in the past to set a robust 



lower limit on the total masses of galaxy clusters (11281) . For an ideal gas equation 



of state the first law of thermodynamics can be solved to yield the specific entropy, 

—2/3 

s = (3^B/2/im a )ln(7 , pg as ) + constant. It is conventional for studies of the hot gas 
in galaxies and clusters to remove the logarithm and the numerical factor to define 

—2/3 

the quantity, S = (k% / jim^T p gas , which we will refer to as the "entropy" in our 
discussion. 

If parameterized models are assumed for both the entropy and mass profiles, then 
the equation of hydrostatic equilibrium can be rewritten so that, 

dJ = _2GM^ll S -V^ ^ p2 /5. (12) 

dr 5 r z 
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Solving this equation for | gives the profiles of gas density, p gas = (P/S) 3 / 5 = 
S -3 / 5 ^ 3 / 2 , and temperature, k^T / [lm a = S 3 / 5 P 2 / 5 = S 3 / 5 ^, which are compared 
to the observation to constrain the parameters of the input S and M(< r) models. 
Since, however, the total mass M{< r) also includes the gas mass which depends 
on |, equation dT2l > can be solved for £, by direct integration only if the gas mass 
is negligible with respect to the mass of the stars and DM. To include the gas mass 
self-consistently, equation ([T21 can be differentiated and rearranged to yield, 

d + ^G r y^ n _ JJ^G + (13) 



dr \ dr 

where use has been made of the relation dM(< r)/dr = 47rr 2 (p sta rs + Pdm + Pgas)- 
The two boundary conditions for this second-order differential equation require 
specifying the values of t, and d% /dr at some radius. The value of B, amounts to 
specifying the pressure at some radius, which is just a normalization constant deter- 
mined by the fitting of the density and temperature profiles. As discussed by J90I) . 
because the gas mass makes a negligible contribution to the total mass at small 
radius, equation (TT2l i can be evaluated at such a small radius to give the d^/dr 
boundary condition. 

The Schwarzschild criterion for convective stability requires that the entropy in- 
creases monotonically with increasing radius, consistent with observations of re- 
laxed elliptical galaxies, for which the entropy profiles are well described by simple 
models of a constant with a broken power-law (9(1 91 ; 92). The Schwarzschild cri- 



terion also effectively limits the magnitude of any temperature gradient. Stronger 
restrictions on the tem perature gradient potentially could be enforced by MTI/HBI 
instabilities (EH EH that develop in the weakly magnetized plasma characteristic 



of the hot gas in a galaxy cluster. These instabilities arise when heat and momentum 
transport occurs only along magnetic field lines, which may happen in regions of 
a plasma where the gyroradius is much smaller than the the Coulomb mean free 
path (such as for an elliptical galaxy - see ^2. It . But since these instabilities are 
expected to be weaker for the lower temperatures o f elliptical galaxies, and they 



are also suppressed by small amounts of turbulence (11311 11321) . their relevance to 
elliptical galaxies is unclear. 

Whereas gas density and temperature profiles display a large range of behavior, 
a robust result of cosmological simulations consideri ng only gr avity is that the ra- 



dial entropy profile follows a power-law, S(r) ~ r l - 1 d 133b 1 1341) . Since the cooling 
and feedback processes that would modify this relation are expected to be most sig- 
nificant closer to the centers of halos, the S(r) ~ r relation provides a physically 
well-motivated asymptotic relation that is particularly useful when extrapolating the 
hydrostatic model outside the data range (e.g., computing model projections). 
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Fig. 3 Chandra temperature profiles of two galaxies from d9l|) . The top panel shows the centrally 
rising profile of NGC 1332 characteristic of galaxies with M v i r < 10 13 M Q and average temperature 
k%T ~ 0.5 keV. The bottom panel shows the centrally falling profile of NGC 4261 characteristic of 
galaxies withM v j r > 10 13 M. ; and average temperature k%T ~ 1 keV. The A'-band half-light radii 
are 2.7 and 3.4 kpc, respectively, for NGC 1332 and NGC 4261. 



3 Is DM Required? 



Prior to the launches of Chandra and XMM, the principal limiting factor for most 
X-ray studies of DM in elliptical galaxies was the temperature profile. (Non- 
spherical DM constraints, which are largely insensitive to the tempera ture profile, 
can be more affected by unresolved discrete sources (11381 ll 39c 1 1401) 



see 



In most cases isothermality was assumed, but the possibility of large radial tem- 
perature gradients meant that no firm conclusions on DM could be obtained for 
most galaxies (14 81) . For a small number of the brightest, hottest (kgT ~ 1 keV) 
most massive (few x 1O 13 M0) galaxies, ROSAT was able to place interesting con- 
straints on t he tem perature pro files, from which large amou nts of DM were i nferr ed: 
NGC 1339 (fl4ll) . NGC 4472 (EUl EHb, NGC 4649 (fl43l and NGC 5044 (flOol 
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Fig. 4 Representative total mass-to-light (M/L) ratios of elliptical galaxies obtained from Chan- 
dra and XMM observations clearly showing rising M/L, and thus the need for DM, outside of 
1 — 2 optical half-light radii. (Top Left) Clear confirmation of rising M/L in the massive galaxies 
NGC 1399 and NGC 5044 (H) found previously with ROSAT observations. (Top Right) Results for 
the galaxies displayed in Figure [T] obtained by i94h . The black regions correspond to the galaxies 
NGC 720, NGC 4325, and NGC 6482 that have M vir < 10 13 M Q while the red regions correspond 
to the galaxies NGC 1407, NGC 4472, NGC 4649, and NGC 4261 that have M vir > 10 13 M & . 
Each resu lt shown r eflects a la confidence region. (Bottom Left) 1 a error regions for three galax- 
ies from | |135| ; Il36t) . where isothermal temperature profiles have been adopted for NGC 57 and 
NGC 7796. (Bottom Right) Results from jl37i) for M/L profiles of galaxies with temperature pro- 
files that fall toward the center (red lines) and for others (blue lines). 



Now accurate temperature profiles have been measured with Chandra and XMM 
for many galaxies (Figure[3]i. The shapes of the radial temperat ure p rofiles of early- 
type galaxies can be separated into roughly two classes (182tl94tll37l) : (1) those with 
profiles that rise toward the center, corresponding to galaxies with M v i r < 10 13 M Q 
and average temperature k%T ~ 0.5 keV, and (2) those with profiles that fall toward 
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R (kpc) 

Fig. 5 Radial mass profile of NGC 720 obtained by f92h using deep observations with Chandra 
and Suzaku. The solid (black) line indicates the total enclosed mass, the dashed (red) line indicates 
the stellar mass, the dotted (blue) line is the DM contribution, and the dash-dot (magenta) line is the 
gas mass contribution. The grey shaded regions indicate the 1 o~ error on the total mass distribution 
obtained using the entropy-based forward-fitting method. Overlaid are a set of data points derived 
from an independent analysis using the traditional smoothed inversion approach; i.e., the model 
indicated by the shaded region was not obtained by fitting these data points. Notice that DM is 
clearly required for radii larger than ~ 5 kpc, about 2R e , where R e = 3.1 kpc is the K-band half- 
light radius. 



the center, corresponding to galaxies with M v i r <; 1O 13 M and average temperature 
k#T ~ 1 keV. These temperature profiles enable accurate DM constraints for many 
galaxies, including lower mass systems (few x 1O 12 M0). 

Equipped with accurate temperature profiles from Chandra and XMM, precise 
con straints on mass-t o-light ratios (M/L ) have n ow b een obtained for many galax- 
ies {HI [HI [H HI I65t ll3d: QH [l4l S SI ElS 0); here M is the total mass 
and L is the stellar luminosity, each computed within the same radius. Generally, for 
those galaxies with good constraints on the temperature profile out to a large radius, 
the data exclude a constant M/L profile at high significance, instead requiring that 
M/L rises with increasing radius outside of 1-2 optical half-light radii (Figure |4j. 
The failure of the constant M/L model provides strong evidence for DM. 

It should be stressed that the current generation of X-ray satellites is able for the 
first time to provide excellent constraints on DM in lower mass galaxies with aver- 
age temperature kWT ~ 0.5 ke V. Combining data from deep observations with Chan- 
dra and Suzaku, d92l) obtained strong constraints on the mass profile of NGC 720 
(see Figure|5]l which was found to have a total mass (M2500 = 1 .6 ± 0.2 x 1O 12 M ) 
similar to that of the Milky Way. 
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This X-ray evidence for DM determined individually for a relatively small 
number of elliptical galaxies (~ 15, depending on the upper mass limit) is com- 
ple mented by evidence from large statistical studies in the o ptical . Weak lens- 
ing ( 147 : 148 : 1491 : 150l) and the dynamics of satellite galaxies ( 151) consider the 
accumulated data from a very large number of galaxies and find that the average 
halo mass ~ 10 12 Mq clearly exceeds that which can be attributed to the stars. As 
for constraints on individual systems, strong evidence for DM using sophisti cated , 
orbit-based stellar dynamical models has been found recently for NGC 4649 (152), 
although the mass somewhat exceeds that inferred by X-rays, a point to which we 
shall return in more detail in SJS] 



4 Radial DM Profile 

Although the need for DM is now established for many elliptical galaxies, the X-ray 
data do not as of yet clearly favor a specific model for the radial profile. In particular, 
both the NFW and ~ r~ 2 DM density profiles are consistent with available data. 
The NFW profile provides a useful means to compare observed DM halos to those 
predicted by cosmological models. The two parameters which specify the NFW 
profile are the concentration and halo mass . The concentration is defined as, 
C A = r A/ r s, where r s is a scale radius (denoting the radius at which the logarithmic 
denstiy slope is -2), and r& is a reference radius defined so that the average density 
of the halo within the sphere of radius r& equals the number A times the critical 
density of the Universe. Quoted values for A typically range from 100-2500. The 
mass, M4, is the mass enclosed within r& . Often A is chose n to b e the solution to the 



spherical collapse model at the time of cluster virialization (11531) . where A ~ 100 at 
the present epoch in the standard ACDM cosmology, in which case the parameters 
are sometimes referred to as the "virial radius", r v i r , "virial mass", M v i r , and "virial 
concentration", c v „. 

Dissipationless simulations of the ACDM model predict c vn m 10 for DM ha- 
los with an intrinsic scatter of sa 0.1 d ex over the mass range of elliptical galax- 



ies considered in this review (I155t 11561) . Some early measurements obtained much 



larger concentrations on these mass scales inconsistent with theory. Firstly, (11571) 
attempted to deconvolve the large, energy-dependent PSF of ASCA observations, 
which is challenging given that the size of the PSF (> 1') is generally as large or 
larger than the optica l half-light radii of nearby elliptical galaxies observed in X- 



rays. Secondly, ( 11581) used ROSAT observations, but they had to assume the gas is 



isothermal t hrou ghout each galaxy. Finally, one early study using Chandra data of 



NGC 6482 (1159b with an accurately measured (non-isothermal) temperature profile 



also obtained a very large concentration (C200 = 61,c v i r ~ 80). 



Mamon & Lokas d 1 541) proposed that such large concentrations disagreeing with 
theoretical expectation were primarily an artifact of neglecting to include a com- 
ponent for the stellar mass in X-ray determinations of the DM profile. This bias 



was confirmed with Chandra observations by (1941) in their sample of seven galaxies 
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Fig. 6 Cyi! — M v ir results (lcr) for the seven elliptical galaxies shown in Figure [T]from (94) where 
the mass model includes (Left Panel) only a single NFW component and (Right Panel) an NFW 
model for the DM and a Hernquist model for for the stellar mass. The omission of the stellar 
co mpon ent leads to a large overestimate of the concentration parameter as originally suggested 
by jl54l) . When including components for both the stars and DM the values are mostly consistent 
with the theoretical expectation within the errors (solid line is median and dashed lines Iff intrinsic 
scatter), with a weak suggestion that the galaxies lie systematically above the theoretical relation. 



(Figure [6), demonstrating that reliable X-ray measurements of the DM concentra- 
tion on the galaxy scale definitely require the stellar mass to be modeled accurately. 
There is good agreement in the concentratio ns ob tained for studies of galaxies like 
NGC 1407 that include the stellar mass (194c 1 1461) . But statistically significant sys- 
tematic differences can arise in some (higher mass) galaxies between studies that 
include stellar mass (165 ) and those that instead exclude a large portion of the cen- 
tral region (66). As emphasized by (65), since an accurate determination of the DM 
concentration requires that the NFW scale radius be accurately measured, it is es- 
sential that r s lie within the range of data being fitted. One the other hand, excluding 
the central region of a galaxy disturbed by AGN feedback is a reasonable approach 
for hydrostatic analysis. In either case, the systematic error in the adopted analysis 
choice needs to be investigated and quantified. 

The observed c vlr — M v - lr relation for elliptical galaxies inferred from X-rays (Fig- 
ure [6} is broadly consistent with the cosmological prediction within the estimated 
observational errors. Despite the relatively large error regions, these measurements 
on the elliptical galaxy scale, when combined with measurements for groups and 
massive clusters, were crucial to providing a sufficiently wide mass baseline to es- 
tablish that the c v h- — M v h- slope is indeed negative (a robust predictio n of C DM mod- 
els) with a value (—0.17 ±0.03) close to the theoretical prediction (11601) . However, 
even when accounting for the stellar mass the normalization of the c v ; r — M v ; r rela- 
tion on the galaxy scale remains larger than predicted. While the effect is not hi ghly 
significant for the results shown in Figure [6] the recent analysis of NGC 720 ( 192b 
using deeper Chandra data mentioned previously with an improved entropy-based 
analysis finds c v h- = 26 ±5, M v h- = 3.3 ±0.4 x 10 12 M Q , which exceeds the theo- 
retical prediction (considering intrinsic scatter) by about 2a. The systematic offset 
could represent one or more of the following: the selection of preferentially early 
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Fig. 7 Gravitational potential profiles and power-law fits from (Il62l) . Chandra data are shown in 
red and XMM in blue. Solid vertical lines indicate the radial range included for the power-law fits, 
while the vertical dotted lines indicate the optical half-light radius. 



forming systems (NGC 720 is a very relaxed isolated galaxy), departures from hy- 
drostatic equilibrium (underestimate of stellar mass contr ibuti on), or the need for a 
different value o f (7r or w in the cosmological model ("see Il60l) . Although adiabatic 



contraction (12.7b 11611) can lower the inferred concentrations, the estimated modest 



reductions still exceed the theoretical prediction (I94t l65b l92l). 

While the NFW profile is a good description of the DM profiles in elliptical 
galaxies, evidence continues to accumulate that the total gravitating mass is very 
well approximated by a single power-law over a wide radial range. Prior to Chan- 
dra and XMM ma ny X-ray studies found mass profiles consistent with p ~ r~ a , 
where a « 2 (see (1162c 1 1631) . and references therein), although isothermal gas was 
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Fig. 8 Radial mass profiles and assoc iated power-law fits. (Left Panel) Results for 10 objects 
spanning galaxies to clusters from ( 163), arbitrarily scaled for clarity. The solid lines are the best- 
fitting profiles determined from an entropy-based forward-fitting procedure, while the data points 
are determined from the traditional smoothed inversion approach; i.e., the displayed models are 
not fi tted t o the data points but are derived i ndep endently. (Right Panel) Results for 6 objects 
from based on X-ray data analyzed by ( 162). Solid lines show the mass best-fitting mass 

profiles, while the dotted lines show the best-fitting power-laws to these mass profiles. 




Fig. 9 From dl63l> : (Left Panel) X-ray (black circles) measurements of logarithmic mass slope 
plotted against optical half-light radius. The red lines are model predictions decomposing a power- 
law total mass profile into Sersic stellar light and NFW DM profiles (see text). (Right Panel) Ratio 
y m of total mass to stellar mass enclosed within R e plotted vs. R e . The solid line is the model 
prediction corresponding to the n = 4 curve shown in the left panel which agre es well with the data 
points corresponding to the ^f-band M/L ratios adapted from the data of | |164|) . 



generally assumed. These results have now been confirmed and strengthened (Fig- 
ures |8][wifl]^c£urate temperature profiles measured with Chandra and XMM 
( I82tl94tll62ill63l : ll 181) . Furthermore, using Chandra data of 10 galaxies span ning a 
wide range in M v j r (including some clusters) and optical half-light radii (R e ), d 163b 
discovered that the power-law index a decreases systematically with increasing R e 
(Figure|9]l. This behavior can be explained by the combination of the stellar (Sersic) 
and DM (NFW) profiles required to produce a power-law total mass profile, which 
leads to the relation a = 2.31 — 0.541og 10 /? e for a Sersic index n w 4 corresponding 
to a stellar de Vaucouleurs profile. 
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This correlation has recently been confirmed by (11651 ) using combined central 
velocity dispersions and strong gravitational lensing in 73 early-type galaxies. These 
optical results provide very important supporting evidence, but since they essentially 
determine a mass profile from two mass data points (i.e., at the center and at the 
Einstein radius of a galaxy), the X-ray measurements remain of cr itical importance. 
This relation between a with R e has far-reaching implications, as (11631) showed that 
it implies a DM fraction within R e that varies systematically with the properties of 
the galaxy in such a manner as to reprod uce, without fine tuning, the observed tilt 
of the Fundamental Plane. Consequently, d 1631) speculated that establishing a nearly 
power-law total mass distribution is a fundamental feature of galaxy formation and 
the primary factor which determines the tilt of the fundamental plane. 



5 Baryon Fraction 

Please see the full article to be published in the book, Hot Interstellar Matter in 
Elliptical Galaxies, eds. D.-W. Kim & S. Pellegrini, Astrophysics & Space Science 
Library (ASSL), Springer. 



6 Non- Spherical Constraints 

Please see the full article to be published in the book, Hot Interstellar Matter in 
Elliptical Galaxies, eds. D.-W. Kim & S. Pellegrini, Astrophysics & Space Science 
Library (ASSL), Springer. 



7 SMBHs 

Please see the full article to be published in the book, Hot Interstellar Matter in 
Elliptical Galaxies, eds. D.-W. Kim & S. Pellegrini, Astrophysics & Space Science 
Library (ASSL), Springer. 



8 Accuracy of the Hydrostatic Equilibrium Approximation 

Provided other sources of systematic error are controlled, the accuracy of a mass 
measurement from X-rays reflects the degree to which non-thermal motions are 
present in the hot gas. As remarked in Sj2T|only limited direct constraints on gas 
motions are possible with current X-ray detectors, and such measurements that exist 
suggest subsonic non-thermal motions. An indirect assessment of non-thermal mo- 
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Fig. 10 Comparison of the stellar M/L ratio s from X -ray mass model ing w ith those predicted 
from three sets of SSP models (Squares: d 167l> ; circles: dl68l) ; triangles: | |169|) ). The dotted line is 
"y=x". The good agreement implies deviations from hydrostatic equilibrium are less than ~10- 
20%. Adapted from |9l|). 



tions can be achieved by comparing the X-ray-determined stellar M/L ratios with 
those predicted for each galaxy by single-burst stellar population synthesis (SSP) 
models. For a sample of 4 galaxies (including NGC 4472 and NGC 4261, both of 
which have AGN-blown cavities in the ISM), the X-rays and SSP display good over- 
all agreement (Fig [Tot, implying non-thermal pres sure is no larger than -~20% (1911) . 
Similar agreement was found for NGC 1407 (11461) . While the SSP models could be 
in error due to a mixture of differently aged st ars in each galaxy or the initial mass 
function (IMF) differing from that of Kroupa (11661) . for non-thermal pressure to be 
significant, it must exist in a finely balanced conspiracy with the shape of the IMF. 

Alternatively, if the three-dimensional gravitational potential can be measured 
independently, it should be possible to use the measured X-ray properties of the gas 
to make an (indirect) constraint on non-thermal gas motions. Since elliptical galax- 
ies with X-ray data of sufficient quality for detailed mass analysis are necessarily 
quite nearby (within ~ 100 Mpc), the only techniques (aside from X-ray tools) 
currently available to measure the gravitational potential with sufficient accuracy 
are axisymmetric (or triaxial) stellar dynamical methods. The most sophisticated of 
these methods employ the full information in the line-of-sight stellar velocity distri- 
bution to break the mass-orbit degeneracy ( 1 70t ll 7 1 J 172 ; 173 ), while making mini- 
mal assumptions about the orbital structure ( 174[" 175b . The fundamentally ill-posed 
nature of the inversion problem at the heart of these "orbit-based " me t hods can be 
mitigated by regularization or a maximum entropy constraint (e.g. [175; fl76l[l7l . 

Only a few studies have attempted to combine the X-ray and optical methods in 
this way. Assuming spherical symmetry and severely limiting the orbital structure 
by employing "Jeans modeling", a modest disagreement between the m ass inferre d 



from the two techniques has been reported for the galaxy NGC 4472 (1177b 11781) 
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While this may indicate on average ^20-50% non-thermal pressure support, it 
could also refle ct the limitations of Jeans modeling to reco ver th e orbital structure 
accurately (e.g.. Il79h . In contrast, using a similar approach d 1 80b argued the X-ray 
method overestimates the mass distribution in the elliptical galaxy NGC 1407 by 
^70%, which is difficult to understand unless the gas is globally out-flowing. The 
exact magnitude of this discrepancy, however, is very uncertain due to large system- 
atic errors associated with assumptions in both the X-ray and the optical analysis. 
(A similar conclusion obtained for the galaxy NGC 3379 (83) is also very uncer- 
tain since unresolved point sources dominate the X-ray emission in that particu- 
lar galaxy d84t 1851) .) Large systematic uncertai nties have also characterized similar 
studies of 6 galaxies undertaken by d88l) and (162), who concluded that, on aver- 
age, the X-ray data are consistent with modest, but highly uncertain, non-thermal 
support (see Figure [XT}. They conclude that, although they find ~30% non-thermal 
pressure in the systems they studied, "f/ze uncertainties in [their] model assumptions 
( e.g., spherical symmetry) are sufficiently large that the contribution [of non-thermal 
pressure] could be consistent with zero." d88l) . By an extension of this argument, it 
could also be significantly larger. 
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r (arcsec) 



Fig. 11 (Left Panel) Comparison of circular velocities for six nearby galaxies inferred from hy- 
drostatic X-ray models (horizontal axis) and from a suite of simple stellar dynamical estimates 
(vertical axis). The observation that the galaxies tend to lie above the dotted line representing equal 
X-ray and optical values implies that the X-ray metho d ten ds to underestimate the mass by around 
30%, implying significant non-th erma l gas motions | |162|) . (Right Panel) Stellar dynamics/X-ray 
comparison for NGC 4649 from 1 11521) . The black lines represent the total mass profile (and Iff 
range) from the stellar dynamics. The green line is a representation of the mass profile inferred 
from X-rays i9(t) . and the red line is the stellar mass. 



A major roadblock to the use of more rigorous axisymmetric (or triaxial) "or- 
bit based" methods for modeling the optical data in such work has been the pro- 
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hibitive computational cost of self-consistently incorporating a DM halo i n suc h a 
scheme. R ecent ly, however, such models have begun to be developed (e.g.. ll8ll) . In 
particular, (11521) found a discrepancy between the mass profile of NGC 4649 from 
orbit-based models and the (spherical) X-ray mass distribution, implying on aver- 
age ~40% non-thermal support. This is significantly larger than was required by 
the hydrodynamical models of (14 II) to reconcile the stellar and X-ray ellipticity pro- 
files (although this comparison is not strictly self-consistent as the hydrodynamical 
mode ls ass umed a gravitational potential consistent with the X-ray results). Simi- 
larly, d 1 731) fou nd d iscrepancies for the galaxy M 87, while the X-ray mass profiles 
were found by (II 18b to underestimate by ^20-45% the gravitating mass determined 
from orbit-based stellar models in the centers of three galaxies, with the sign of the 
discrepancy reversing at large (> 10 kpc) radii. While inclination effects could con- 
ceivably contribute a significant uncertainty into the stellar dynamical analysis of 
the round galaxies considered in these studies, a similar mismatch between the X- 
ray and (orbit-based) optical mass estimates (corresponding to ^50% non-thermal 
support) ha s bee n reported in the central ^3 kpc region for the edge-on SO galaxy 
NGC 1332 (EH. 

Hence, the current status of comparisons between stellar dynamics and X-rays is 
that X-rays tend, on average, to underestimate the mass obtained by stellar dynam- 
ics by ~ 30%, but with considerable uncertainty on many of the measurements, and 
the comparison can be complicated by a radial dependence in the discrepancy. To 
better understand the systematic errors in these comparisons, it is essential in the 
future that self-consistent comparisons be made; i.e., simultaneous modeling of the 
X-ray and stellar dynamical data. The X-ray data analysis is subject to numerous 
potential systematic effects which need to be controlled; e.g., background charac- 
terization, measurements of the metal abundances, treatment of unresolved discrete 
sources, and even unce rtainties associated with the mass modeling technique that is 



adopted (e.g. . I91tll 181) . While many of these systematic errors are modest in abso- 



lute terms, they can be comparable to, or at times exceed, the statistical errors. With 
the detailed comparison between X-ray and optical mass profiles now demanding 
the interpretation of differences as small as ^20%, precise control of these errors 
is essential. Until this becomes routine, the accurate radial profiles of possible non- 
thermal gas motions inferred by comparing mass profiles determined from X-rays 
and stellar dynamics will remain out of reach. 



9 Conclusions 

With the advent of accurate radial temperature profiles from Chandra, XMM, and 
Suzaku, the need for DM is now firmly established in many giant elliptical galaxies 
(1O 12 M < M vir < 10 13 M o ). Current data do not distinguish between NFW and 
isothermal DM profiles. For the cosmologically motivated NFW model the inferred 
concentration parameters of elliptical galaxies generally exceed the mean theoreti- 
cal relation. This discrepancy may be due in large part to the preferential selection 
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of the most relaxed galaxies for X-ray studies. X-ray observations confirm that the 
total mass profile (baryons+DM) is close to isothermal (M ~ r), and new evidence 
suggests a more general power-law relation for the slope of the total mass profile 
that varies with the stellar half-light radius. While the global baryon fraction re- 
mains difficult to constrain, a recent joint Chandra-Suzaku analysis of NGC 720 
indicates a global baryon fraction consistent with the cosmological value, suggest- 
ing a Milky Way-size galaxy can retain all of its baryons. The axial ratio of the DM 
inferred from the X-ray isophote shapes of NGC 720 is very consistent with those of 
halos produced in cosmological simulations. Finally, the unprecedented high spatial 
resolution of Chandra has enabled the first measurements of the mass of SMBHs 
in a few elliptical galaxies, obtaining a precision comparable to that achieved by 
traditional optical methods in NGC 4649. 

The most important uncertainty in X-ray determinations of DM in elliptical 
galaxies is no longer the measurement of the temperature profile but now the ac- 
curacy of the hydrostatic equilibrium approximation; i.e., the magnitude of non- 
thermal motions in the hot ISM. The limited direct constraints possible with current 
data suggest subsonic non-thermal motions, and improved measurements of turbu- 
lent motions via Doppler line broadening in the brightest systems are expected from 
Astro-H in a few years. In the meantime, it is essential to perform self-consistent 
studies of X-ray and stellar dynamical data of elliptical galaxies to assess the mag- 
nitude of non-thermal gas motions and achieve more reliable DM measurements 
than possible with either technique individually. 
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